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ABSTRACT 

We investigate the long-term dynamical stability of hypothetical moons orbiting extrasolar giant planets. 
Stellar tides brake a planet's rotation and, together with tidal migration, act to remove satellites; this process 
limits the lifetimes of larger moons in extrasolar planetary systems. Because more massive satellites are removed 
more quickly than less massive ones, we are able to derive an upper mass limit for those satellites that might 
have survived to the present day. For example, we estimate that no primordial satellites with masses greater than 
7 x 10 -7 M ffi (~70 km radius for p = 3 g cm~ 3 ) could have survived around the transiting planet HD209458b 
for the age of the system. No meaningful mass limits can be placed on moons orbiting Jovian planets more than 
~ 0.6 AU from their parent stars. Earth-like moons of Jovian planets could exist for 5 Gyr in systems where the 
stellar mass is greater than 0.15 M Q . Transits show the most promise for the discovery of extrasolar moons — 
we discuss prospects for satellite detection via transits using space-based photometric surveys and the limits on 
the planetary tidal dissipation factor Q p that a discovery would imply. 
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Each of the giant planets in our solar system posesses a 
satellite system. Since the discovery of planets in other so- 
lar systems (Marcy et al. 2000b) , the question of whether 
these extrasolar planets also have satellites has become 
relevant and addressible. Extrasolar planets cannot be 
observed directly with current technology, and observing 
moons around them poses an even greater technical chal- 
lenge. However, high precision photometry of stars during 
transits of planets can detect extrasolar moons either by di- 
rect satellite transit or through perturbations in the timing 
of the planet transit (Sartoretti & Schneider 1999). Using 
these techniques, Brown et al. (2001) placed upper limits 
of 1.2 Earth radii ( R®) and 3 Earth masses ( M®) on any 
satellites orbiting the transiting planet HD209458b based 
on the Hubble Space Telescope transit lightcurve. 

The tidal bulge that a satellite induces on its parent 
planet perturbs the satellite's orbit (e.g., Burns 1986), 
causing migrations in scmimajor axis that can lead to the 
loss of the satellite. For an isolated planet, satellite re- 
moval occurs either through increase in the satellite's or- 
bital semimajor axis until it escapes, or by inward spiral 
until it impacts the planet's surface (Counselman 1973). In 
the presence of the parent star, stellar-induced tidal fric- 
tion slows the planet's rotation, and the resulting planet- 
satellite tides cause the satellite to spiral inward towards 
the planet (Ward & Reid 1973; Burns 1973). This effect 
is especially important for a planet in close proximity to 
its star, and has been suggested as the reason for the lack 
of satellites around Mercury (Ward & Reid 1973; Burns 

1 07.^ 



In this paper, we apply tidal theory and the results of 
numerical orbital integrations to the issue of satellites or- 
biting close-in extrasolar giant planets. We place limits on 
the masses of satellites that extrasolar planets may posess, 
discuss the implications these limits have for the detection 
of extrasolar satellites, and apply our results to the issue 
of Earth-like satellites orbiting extrasolar giant planets. 

2. TIDAL THEORY AND METHODS 

According to conventional tidal theory, the relative val- 
ues of the planetary rotation rate, fl p , and the orbital mean 
motion of the moon, n m (both in units of radians/sec), de- 
termine the direction of orbital evolution (see e.g., Murray 
& Dermott 2000). For a moon orbiting a planet slower 
than the planet rotates (n m < fip), the tidal bulge in- 
duced on the planet by the satellite will be dragged ahead 
of the satellite by an angle S, with tan(2J) = 1/Q p . Here, 
Q p is the parameter describing tidal dissipation within the 
planet (after Goldreich & Soter 1966), with 1/Q P equal to 
the fraction of tidal energy dissipated during each tidal cy- 
cle. Gravitational interactions between the tidal bulge and 
the satellite induce torques that transfer angular momen- 
tum and dissipate energy, slowing the planet's rotation and 
increasing the orbital semimajor axis of the satellite. Con- 
versely, for satellites orbiting faster than their planet's ro- 
tation (n m > Hip), the planet is spun up and the satellite's 
semimajor axis decreases. The same mechanism causes 
torques on the planet from its parent star which slow the 
planet's rotation (Murray & Dermott 2000). 

The torque on the planet due to the tidal bulge raised 



by the moon (r p _ m ) is (Murray & Dermott 2000) 
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where k 2p is the tidal Love number of the planet, R p is the 
radius of the planet, and G is the gravitational constant. 
The term sign(fi — n m ) is equal to 1 if (O — n m ) is pos- 
itive, and is equal to —1 if it is negative. We obtain the 
expression for the stellar torque on the planet by replacing 
M m , the mass of the moon, with M*, the mass of the star; 
by replacing a m , the semimajor axis of the satellite's or- 
bit with dp, the semimajor axis of the planet's orbit about 
the star (circular orbits are assumed); and by using the 
planet's mean orbital motion n p instead of n m : 
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The moon's semimajor axis, a m , and the moon's mean 
motion, n m , are related by Kepler's law, n^Bp = GM p . 
These torques affect both n m and Q p . The rate of change 
of ftp is obtained by dividing the total torque on the planet 
by the planet's moment of inertia: 



dfl p 
~dt 
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where I p is the planet's moment of inertia. 



Under the circumstances studied in this paper, where a 
planet is orbited by a much smaller satellite, t p _* is much 
greater than r p _ TO for most of the system's lifetime. Be- 
cause the moon's orbital moment of inertia depends on 
n m , the equivalent expression for n m is less trivial to de- 
rive. We obtain it by setting the torque equal to the rate 
of change of the angular momentum and solving for dn/dt 
using the planet's mass, M p (e.g., Peale 1988): 
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Given appropriate initial and boundary conditions, inte- 
gration of Equations 3 and 4 determines the state of the 
system at any given time. 

An important boundary condition for such an integra- 
tion is the critical semimajor axis, or the location of the 
outermost satellite orbit that remains bound to the planet. 
This location must be within the planet's gravitational in- 
fluence, or Hill sphere, and has been generally thought to 
lie between 1/3 and 1/2 the radius of the Hill sphere (Rh) 
(Burns 1986), where 
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Recently, Holman & Wiegert (1999) investigated the sta- 
bility of planets in binary star systems and their results 
are applicable to the planet-satellite situation as well. 
Through numerical integrations of a test particle orbiting 



(1999) found that for high mass ratio binaries, the criti- 
cal semimajor axis for objects orbiting the secondary in 
its orbital plane is equal to a constant fraction (/) of the 
secondary's Hill radius, or 
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We treat a star orbited by a much less massive planet as a 
high mass ratio binary system and deduce that the critical 
semimajor axis for a satellite orbiting the planet is .36Rh 
(f = .36) for prograde satellites (from Holman & Wiegert 
1999, Figure 1). This agrees closely with Burns (1986). 
In fact, none of the prograde moons of our solar system 
orbit outside this radius (see Table 1). Holman & Wiegert 
(1999) did not treat objects in retrograde orbits (which 
are expected to be more stable than prograde ones), so to 
treat possible captured satellites we take /retrograde = 0.50 
based on the solar system values for a m /RH in Table 1. 

3. CONSTRAINTS ON SATELLITE MASSES 

Satellites orbiting close-in giant planets fall into one 
of three categories based on the history of their orbital 
evolution. Satellites that either start inside the planet's 
synchronous radius (the distance from the planet where 
n m = Qp) or become subsynchronous early in their life- 
times, as a result of the slowing planetary rotation, spend 
their lives spiraling inward toward the cloud tops. Eventu- 
ally these moons collide with the planet or are broken up 
once they migrate inside the Roche limit. Moons that start 
and remain exterior to the synchronous radius evolve out- 
ward over the course of their lives and, given enough time, 
would be lost to interplanetary space as a result of orbtial 
instability. In between these two is a third class of orbital 
history. In this case, a satellite starts well outside the 
synchronous radius and initially spirals outward, but its 
migration direction is reversed when the planet's rotation 
slows enough to move the synchronous radius outside the 
moon's orbit. These moons eventually impact the planet. 

In order to determine which satellites might still exist 
around any given planet, we determine the maximum life- 
time for a moon with a given mass in each orbital evolution 
category. Inward-evolving satellites should maximize their 
lifetime by starting as far from the planet as possible, at 
the critical semimajor axis a CT it (Equation 6), and spiral- 
ing inward all the way to the planet. Outward-evolving 
satellites can survive the longest if they start just outside 
the synchronous radius of the planet, then spiral outward 
to the critical semimajor axis. The maximum lifetime for 
the out-then-in case occurs when a satellite reverses migra- 
tion direction at the outermost possible point, the critical 
semimajor axis. In this case, the moon starts at the semi- 
major axis that allows for it to have reached a cr u by the 
time its planet's synchronous radius also reaches a cr it, thus 
maximizing the time for its inward spiral (see Figure 1). 
For a given satellite mass, if the maximum possible lifetime 
is shorter than the age of the system, then such a satellite 
could not have survived to the present. Because the orbits 



be placed on the masses of satellites that could still exist 
around any given planet. 



3.1. Analytical Treatment 

For a given semimajor axis of a moon, a m , the migration 
rate is the same whether the moon is moving inward or out- 
ward (with the assumption that Q p is independent of the 
tidal forcing frequency Q p — n m ), and the migration rate is 
much faster for satellites close to their parent planets. For 
both the inward- and outward-migrating categories, the 
total lifetime of a satellite (T) is well-characterized by the 
time necessary for a satellite orbit to traverse the entire 
region between the critical semimajor axis (a m — a cr ,t) 
and the planet's surface (a m — R p ), (Murray & Dermott 
2000): 
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Since R p <C a C rit and the exponents are large, the R p 
term inside the parenthesis can be neglected. Substitut- 
ing a cr it = JRh (Equation 6), allowing T to be equal to 
the age of the system, and solving for M m collectively re- 
sult in an analytical expression for the maximum possible 
extant satellite mass in both the inward and outward cases, 
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which is the equation of the bottom, dot-dashed line in 
Figure 2. In the case of satellites that evolve outward then 
inward, the spindown of the planet is important. These 
moons can be saved temporarily by the reversal of their 
orbital migration. This reversal prolongs their lifetimes, 
but by less than a factor of two, because the satellite trans- 
verses the region where a m < a cr u twice. The upper mass 
limit for these satellites is 
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and this limit is plotted as the upper dotted line in Figure 
2. 

We obtain the boundaries between the in, out-then-in, 
and out cases by comparing the time necessary to despin 
the planet to the age of the system. The time necessary 
to spin down the planet to the point that the synchronous 
radius becomes exterior to the critical semimajor axis is 
equal to (Guillot et al. 1996) 
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where f2 p o is the initial planetary rotation rate and Q p i is 
the rotation rate at which the planet's synchronous radius 
is coincident with & C rit' 

At this point, from Kepler's law 

we infer that 
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Fig. 1. — Satellite orbital semimajor axis versus time for 
the maximum-mass moon in three different hypothetical 
4.6 Gyr old 1 M Q , 1 Mj up planetary systems. The solid 
lines represent a system with a planet-star distance of 
0.15 AU, with the thick and thin lines corresponding to 
the satellite semimajor axis and planet synchronous radius 
respectively. Tides between this planet and its star spin 
down the planet in short order, and the moon spends the 
majority of its lifetime evolving inward through tidal inter- 
actions with the planet. It is destroyed upon reaching the 
cloud tops of the planet. The dotted lines correspond to a 
system in which the planet-star separation is 0.20 AU; this 
planet is despun in just under half the age of the system. 
Thus its maximum-mass moon initially moves outward due 
to tidal influences, but later reverses direction due to the 
spindown of the planet and eventually crashes into it. This 
moon reverses direction at the critical semimajor axis (the 
outermost stable orbit point) because doing so maximizes 
its orbital lifetime. The dashed lines are for a planet or- 
bital semimajor axis of 0.25 AU. The star's tidal torques 
on the planet have less influence at this distance, and the 
planet does not despin sufficiently over its lifetime to re- 
verse the orbital migration of its maximum-mass satellite. 
This satellite is lost into interplanetary space due to orbital 
instabilities. 
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For the case where T sp i nc t own -C T (where T is the sys- 
tem liftime) , the maximum-mass moon evolves inward and 
Equation 8 should be used. When the system age is greater 
than this spindown time {T spin down ^> T), moons evolve 
outward and again Equation 8 is valid. However, when 
T S pindown ~ T, the reversal of satellite orbital migration is 
important and Equation 9 provides a more robust upper 
mass limit for surviving satellites. 

These results for the maximum M m are limited by the 
requirement that the rate of angular momentum transfer 
between the planet and the satellite must be less than that 
between the planet and the star when a m is relatively large 
(i.e., t p _* > Tp- m in Equation 3) such that synchroniza- 
tion between the planet and moon does not occur. In 
the case of rocky satellites orbiting gaseous planets, this 
condition is met. For large moon-planet mass ratios or 
large a p , this assumption breaks down, yielding a situa- 
tion more closely resembling the isolated planet-satellite 
systems treated in Counselman (1973). In this case the 
planet and moon can become locked into a 1:1 spin-orbit 
resonance with each another, halting the satellite's orbital 
migration and extending its lifetime. For extrasolar Jovian 
planets (0.3M Q < M p < 13.0M Q ) this occurs when satel- 
lite masses become very large, i.e. greater than 8M0 for 
a lMj U p planet. Such a moon is large enough to accrete 
hydrogen gas during its formation, however, and in such a 
case the system is better treated as a binary planet, taking 
into account the tidal torques of each body on the mutual 
orbit. We do not address that situation here. 

We assume prograde, primordial satellites, but objects 
captured into orbit by a planet late in its life could also 
remain in orbit. We do not treat the physics of satellite 
capture, but the lifetimes of such moons would be affected 
by the same processes described above if prograde, and 
limited by inward migration like Neptune's moon Triton 
(McCord 1966) if retrograde. In the case of retrograde, 
captured moons, the following upper limit on their survival 
lifetime can be placed by rearranging Equation 8: 
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Our analysis assumes a single satellite system. Inward- 
migrating moons could not be slowed significantly by en- 
tering into a resonance with another satellite further in 
because the interior satellite would be migrating faster 
(due to the a~ 6 dependance of the torque in Equation 
1), unless its mass is less than .08 M m (assuming a 2:1 
resonance). Slowly- migrating moons exterior to the satel- 
lite in question cannot slow its orbital migration because 
objects in diverging orbits cannot be captured into reso- 
nances. However, outward-migrating satellites could have 
their lifetimes extended by entering into a resonance with 
an exterior neighbor through intcrsatcllitc angular momen- 
tum transfer (Goldreich 1965), similar to the resonances 
currently slowing the outward migration of Io from Jupiter. 
Thus some outward-limited satellites above the mass limit 
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Fig. 2. — Stability diagram for satellites in a hypothetical 
4.6 Gyr old 1 M Q , 1 Mj up planetary system as a function 
of the planet's orbital semimajor axis. The solid line repre- 
sents the results from numerical integrations of Equations 
3 and 4, and the broken lines are the analytical approxima- 
tions given in Equations 8 (lower, dot-dashed) and 9 (up- 
per, dotted). Satellites above this line are excluded, while 
those which lie below the line may or may not still exist 
depending on their specific orbital evolutionary histories. 
For this specific case, below a p <~ 0.15AU satellite masses 
are limited by their inward migration, above a p <~ 0.23AU 
by their outward migration, and in between by outward 
followed by inward migration with the reversal being the 
result of slowing planetary rotation. 
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3.2. Numerical Treatment 

To verify the limits stipulated in Equations 8 and 9, 
we integrate Equations 3 and 4 numerically, from the ini- 
tial rotation rate and semimajor axis until the satellite's 
demise either through impact with the planet or through 
orbital escape. We use an adaptive stepsize Runge-Kutta 
integrator from Press et al. (1992) and have verified that it 
reproduces our analytical results for small satellite masses. 
We also assume that only Q p and a m change over time — 
other planetary parameters such as Q pi a p , R p , and all 
others are taken to be constant for the length of the in- 
tegration. The expected changes in the planet's orbital 
semimajor axis a p over the course of the integration do 
not significantly affect the calculations, and larger plane- 
tary radii R p in the past would only serve to further reduce 
lifetime of a given satellite beyond what we have calculated 
here, pushing the upper surviving satellite mass lower. 

For each planet, we determine the maximum satellite 
mass that could survive for the observed lifetime of the sys- 
tem by optimizing the initial semimajor axis of the satellite 
so as to maximize its lifetime and then tuning the satellite 
mass until this lifetime is equal to the system age. Sam- 
ple evolutionary histories of this maximum mass satellite 
for a hypothetical 1M©, lMj up system from each orbital 
evolutionary history category are shown in Figure 1. The 
numerically determined maximum mass as a function of 
planetary orbital semimajor axis is shown in Figure 2 as 
the solid line. These numerical results are consistent with 
our analytical upper mass limits from Equations 8 and 9. 

4. IMPLICATIONS 

4.1. Known Extrasolar Planets 

In applying these results to the specific test case of 
the transiting planet HD209458b, we adopt the values 
M* = 1.1 M Q , M p = 0.69 M Jup , a p = 0.0468 AU, 
T = 5.0 Gyr (Mazeh et al. 2000), and R p = 1.35 R Jup 
(Brown et al. 2001) based on observational studies. We 
take fc2p for the planet to be 0.51, the value for an n = 1 
polytrope (Hubbard 1984). The least constrained param- 
eter is the tidal dissipation factor Q p ; for HD209458b we 
adopt Q p = 10 5 , which is consistent with estimates for 
Jupiter's Q p (Goldreich & Soter 1966). However, Q p is not 
known precisely even for the planets in our own solar sys- 
tem, and the precise mechanism for the dissipation of tidal 
energy has not been established. Q p for extrasolar planets, 
and especially for ones whose interiors differ from Jupiter's 
such as close-in giant planets (Burrows et al. 2000), may 
differ substantially from this value. 

Because HD209458b was likely tidally spun down to 
synchronous rotation very quickly (Guillot et al. 1996), 
satellites around it are constrained by the infall time and 
we use Equation 8 to obtain an upper limit of 7 x 10~ 7 M© 
for their masses. Assuming a density of p = 3 g/cm 3 , 
the largest possible satellite would be 70 km in radius 



Amalthea. These limits are consistent with the those 
placed on actual satellites observationally by Brown et al. 
(2001). It is possible for captured satellites to exist around 
HD209458b. Their lifetimes, however, would be exceed- 
ingly short — a 1M© satellite could survive for only 
30, OOOyr (Equation 12), making the probability of detect- 
ing one low unless such captures are common. 

We also calculate the maximum masses for surviving 
moons around other detected extrasolar planets; the re- 
sults are in Table 2. We take the planet masses to be 
equal to the minimum mass determined by radial velocity 
monitoring, as the orbital inclination has not been reli- 
ably determined for any planet except HD209458b. We 
use the same ki p and Q p as we did for HD209458b, but 
take R p — — Rjup 

because the radii for these objects is un- 
known. For this table, we have only chosen planets whose 
orbital eccentricities are less than 0.10 because Equation 
6 applies only to planets in circular orbits. The critical 
semimajor axis for planets in noncircular orbits has not 
yet been determined, thus we leave the calculation of up- 
per mass limits for satellites around these eccentric planets 
for future work. 

4.2. Earth-like Moons 

Our approach can also shed light on the issue of Earth- 
like satellites, which we define to be moons capable of sup- 
porting liquid water. Low-mass satellites do not fit this 
definition due to their inability to retain volatiles (Williams 
et al. 1997). Here we note that high-mass satellites may 
not survive for long periods around close-in planets because 
planets with low masses have smaller Hill spheres and 
therefore, for a given satellite mass, also have shorter max- 
imum moon lifetimes. To calculate in general which giant 
planets might harbor Earth- like satellites, we use Equation 
8 and constrain a p based on the insolation at the planet, F, 
relative to the Earth's insolation F© = 1370 W m~ 2 . We 
use the rough approximation (Hansen & Kawaler 1994) 
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for the stellar luminosity, L» , together with the insolation 
at the planet, 
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to fix the planet's semimajor axis by solving for a p . By 
plugging the resulting value for a p into Equation 8, we can 
exclude Earth-like moons around planets in systems that 
don't satisfy the inequality 
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Equation 15 is plotted in Figure 3 for the same values of 
k2 P , Qp, and / as we use for HD209458b, with R p — Rj up 
and T — 5 Gyr. 

AiVilliamc at al (~IQQ7\ fnnnrl tVio lr>wor limit M > 
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Fig. 3. — Stability of Earth-like moons. For stars of a 
given mass, we estimate a luminosity from Equation 13, 
and then calculate the limits on planet masses that allow 
for the stability of moons with insolation F = 1370 W m~ 2 
and mass 1 Mmna (solid line, left), 1 Mmsfs (solid line, 
center), and 1 MEarth (solid line, right). Planets with the 
proper a p (top) and masses above the appropriate solid 
line might still harbor Earth-like satellites after 5 Gyr. For 
this plot, we have taken R p = R,j up , and as such the values 
are not valid for planets with masses less than 0.3 Mj up . 
Planet mass limits for very low mass stars where M p ~ M* 
are also suspect. From this plot we infer that Earth- like 
satellites of Jovian planets are plausibly stable for 5 Gyr 
around most stars with masses greater than 0.15 M©. 



timescales. Using this mass, we find that Earth-like moons 
orbiting Jovian planets could survive for solar system life- 
times around stars with masses greater than 0.15 M©, and 
that Earth-mass satellites which receive similar insolation 
to the Earth are stable around all Jovian planets orbit- 
ing stars with M* > 0.5 M©. Planets with masses less 
than 0.3 Mj up differ in radius, Q p , and interior structure 
from those with masses greater than 0.3 Mj up . In addi- 
tion, for lower planet masses the planet/satellite mass ratio 
increases beyond the assumption of non-synchronization 
between the planet and moon. For these reasons, we do 
not treat the question of Earth-like satellites of ice-giant 
planets (M p < 0.3M Jup ) here. 

The radial-velocity planet most likely to harbor Earth- 
like moons is HD28185b because of its circular, a p = 1AU 
orbit around a star similar to the Sun withspectral type 
G5V and L* = 1.09L© (Santos et al. 2001). Because the 
calculated upper satellite mass for this planet is above 
8 M©, we can not rule out any satellite masses for this 
object. Thus, Earth-like moons with any mass could plau- 
sibly be stable around HD28185b. 

4.3. Future Discoveries 

Several missions to search for extrasolar planet tran- 
sits by high-precision space-based photometry are in the 
planning stages and will, if launched, have the capability 
of detecting satellites (Sartoretti & Schneider 1999). The 
probability that a given planet will transit across its par- 
ent star decreases with planetary orbital semimajor axis as 
1/dp. Hence these surveys will preferentially detect plan- 
ets in orbits close to their parent stars. However, we have 
shown that it is unlikely that these close-in objects will har- 
bor satellites. Therefore satellite transits are most likely to 
be detected around planets orbiting at moderate distances 
from their parent star ( 0.3AU < a p < 2AU), even though 
planet transits are most likely at small orbital distances. 
If a satellite were detected, Equation 9 could be used to 
place limits on the planetary tidal dissipation parameter 
Q p . By using extreme values of the lifetimes, masses, and 
possible values Q p that may exist, we estimate this process 
will not significantly affect planets more than 0.6AU from 
their parent star, leaving any satellite systems they might 
posess intact. 
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Table 1 
Satellite Semimajor Axes 



Planet 


Satellite 


a-m/Rii 


Earth 


Moon 


0.257 


Mars 


Deimos 


0.0216 


Jupiter 


Callisto 


0.0354 


Jupiter 


Elara 


0.221 


Jupiter 


Sinope 


0.446 R 


Saturn 


Titan 


0.0187 


Saturn 


Iapetus 


0.0545 


Saturn 


Phoebe 


0.198 R 


Saturn 


S/2000 S 9 


0.283 


Uranus 


Obcron 


0.00837 


Uranus 


Setebos 


0.352 R 


Neptune 


Triton 


0.003 R 


Neptune 


Nereid 


0.0475 



R Retrograde 

Note. — Orbital semimajor axes 
of selected solar system satellites are 
listed as a function of their their par- 
ent planet's hill sphere radius, Rh- 



Table 2 

Constraints on Satellites around Selected Extrasolar Planets 



1\[ q m o 
IN cLIIIL- 


Star 


1V1 Olll 1 




c 


Max 






iY.L.lL.1 L.I1L-L. 






(Mj ) 

V JW J up } 


CATH 




Moon 




lVTnnn 














Mass 
(M e ) 




xxaui ua 
(Kill) 




HD83443b 


(5) 


0.35 


0.038 


0.08 


8 x 10- 


8 


30 


Mayor ct al. (2000) 


HD46375b 


(5) 


0.25 


0.041 


0.04 


6 x 10- 


8 


30 


Marcy et al. (2000a) 


HD187123b 


(5) 


0.52 


0.042 


0.03 


6 x 10- 


7 


60 


Butler et al. (1998) 


HD209458b 


5 


0.69 


0.045 





7 x 10- 


7 


70 


Fischer et al. (2002) 


HD179949b 


(5) 


0.84 


0.045 





3 x 10- 


6 


110 


Tinney et al. (2001b) 


HD75289b 


(5) 


0.42 


0.046 


0.053 


6 x 10- 


7 


60 


Udry ct al. (2000) 


BD -10 3166 b 


(5) 


0.48 


0.046 


0.05 


7 x 10- 


7 


70 


Butler et al. (2000) 


T Boo b 


2 


4.1 


0.047 


0.051 


5 x 10- 


4 


600 


Butler et al. (1997) 


51 Pegasus b 


(5) 


0.44 


0.051 


0.013 


2 x 10" 


6 


90 


Marcy et al. (1997) 


U And b 


2.6 


0.71 


0.059 


0.034 


2 x 10- 


5 


190 


Butler et al. (1999) 


HD168746b 


(5) 


0.24 


0.066 





2 x 10- 


6 


100 


l 


HD130322b 


(5) 


1 


0.088 


0.044 


0.0008 


730 


Udry et al. (2000) 


55Cnc b 


5 


0.84 


0.11 


0.051 


0.001 




810 


Butler et al. (1997) 


G186b 


(5) 


3.6 


0.11 


0.042 


0.1 




3950 


Queloz et al. (2000) 


HD195019b 


3.2 


3.5 


0.14 


0.03 


0.8 




7090 


Fischer et al. (1999) 


GJ876c 


5 


1.9 


0.21 


0.1 


6 




14050 


Marcy et al. (2001) 


rho CrB b 


10 


1.1 


0.23 


0.028 


0.3 




5310 


Noyes et al. (1997) 


U And c 


2.6 


2.1 


0.83 


0.018 








Butler et al. (1999) 


HD28185b 


(5) 


5.6 


1 


0.06 








Santos et al. (2001) 


HD27442b 


(5) 


1.4 


1.2 


0.02 








Butler et al. (2001) 


HD114783b 


(5) 


1 


1.2 


0.1 








Vogt et al. (2002) 


HD23079b 


(5) 


2.5 


1.5 


0.02 








Tinney et al. (2001a) 


HD4208b 


(5) 


0.8 


1.7 


0.01 








Vogt et al. (2002) 


47UMa b 


6.9 


2.5 


2.1 


0.061 








Fischer et al. (2002) 


47UMa c 


6.9 


0.76 


3.7 


0.1 








Fischer et al. (2002) 



1 http : / /obs www . unige . ch / ~udry / planet /hd 1 68746 . html 
2 http://c3po. lpl.arizona.edu/egpdb 

Note. — Upper satellite mass limits are determined with Q p — 10 5 , R p — Rj up , and with M p = 
msini. Where no system lifetime was available in the literature, we have taken the system age to 
be 5 Gyr, and those cases are indicated by parenthesis. To obtain maximum moon radii, we assume 
p = 3 g cm~ 3 . We can not place useful limits for those planets where maximum masses and radii 
are not listed. Planets with orbital eccentricities greater than 0.1 are excluded due to the difficulty 
in determining the proper value of /. The planet data used to generate this table have been formed 
into a world wide web accessible database of extrasolar planets ? 



